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Abstract 

Causal dynamical triangulations (CDT) constitute a background in- 
dependent, nonperturbative approach to quantum gravity, in which 
the gravitational path integral is approximated by the weighted sum 
over causally well-behaving simplicial manifolds i.e. causal triangu- 
lations. This thesis is an analysis of the data from the Monte Carlo 
computer simulations of CDT in 3 + 1 dimensions. It is confirmed 
here that there exist the semiclassical limit of CDT for so-called (4, 1) 
(or equivalent (1,4)) simplices, being a discrete version of the mini- 
superspace model. Next, the form of the corresponding discrete action 
is investigated. Furthermore, it is demonstrated that the effective, 
semiclassical solution works also after the inclusion of remaining (3, 2) 
and (2, 3) simplices, treated collectively. A specific form of the result- 
ing extended discrete action is examined and a transition from the 
broader framework to the former narrower one is shown. 



*Master Thesis, supervised by prof, dr hab. J. Jurkiewicz and defended in June 2010. 
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Introduction 



In spite of numerous efforts, based on the most diverse assumptions, quanti- 
sation of gravitational field remains an elusive task. Regarding our standard 
approach to quantum field theory, perturbative nonrenormalisability is con- 
sidered the fundamental obstacle here. Therefore attempts to formulate a 
nonperturbative theory of quantum gravity have gained a lot of interest. 
Probably the most intuitive path in this direction is founded on Feynman 
path integral defined for the geometry of spacetime (see Sec. 1.1). 

One of the approaches to quantum gravity with such a starting point 
is known by the name of causal dynamical triangulations (CDT). It is a 
modification (Sec. 1.4) of earlier Euclidean dynamical triangulations (EDT) 
(Sec. 1.2, 1.3) and consists in replacing the gravitational path integral with a 
weighted sum over piecewise linear manifolds assembled from simplices (i.e. 
triangulations), every of which must preserve causality and admit a global 
proper time foliation. Such an approach is background independent, as it 
does not favour any particular geometry. 

CDT can be naturally Wick rotated and becomes then a statistical system 
of random geometry. As a result, it is suited to be investigated numerically 
using Monte Carlo methods, in particular in 3 + 1 dimensions (Sec. 2.1). The 
corresponding program has been written by J. Jurkiewicz and A. Gorlich and 
improved only recently by J. Gizbert-Studnicki, who also set up the latest 
simulations. Hence one should refer to [19] for more details concerning nu- 
merics. The following thesis is an analysis of the obtained data and its point 
of departure is confirmation of the existence of a physically relevant semiclas- 
sical limit of CDT, being the discretised mini-superspace model (Sec. 2.2). 

In subsequent chapters I examine this effective, semiclassical solution of 
CDT employing two kinds of observables: average distributions of the number 
of simplices and the analogical (inverse) covariance matrices. Particularly, 
I seek in the form of the discrete effective action. In Chapter 3 I consider 
only so-called (4,1) simplices, extend the description to (3,2) simplices in 
Chapter 4 and sketch out the possibility of separate inclusion of (3, 2) and 
(2, 3) simplices in Appendix A. The whole analysis has been conducted with 
Mathematica and all graphs have been plotted with gnuplot. 
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Chapter 1 

CDT model of quantum gravity 



1.1 Path integral for gravitation 

Classical field theory of gravitation, known as general relativity is specified by 
the Einstein-Hilbert action, which for a differentiable manifold M of d = D+l 
dimensions has the form 

SehM = Jj d x y/\teti\(R - 2A) (1.1) 

V 87rG JdM J 

where g is a metric 2-form field, R a Ricci scalar, G denotes the gravitational 
constant and A the cosmological constant while the integral in brackets, 
appearing only when M has boundary, is the Gibbons-Hawking- York term 
with an induced metric h and trace of an extrinsic curvature K. 

Formally, one defines Feynman path integral for gravitational field, which 
actually is a grand canonical partition function for quantum gravity, as 

Z{G,A):= [ V[g, u ] eMiS EH [g„A), (1-2) 

JGeom(M) 

where g^ v = [g^ v ] G Geom(M) denotes the equivalence class of metrics with 
respect to diffeomorphisms - a geometry. In particular, the partition function 
for manifold with spatial boundaries at proper times t', t", 

GG,A{9ij,9ij;t",t') := / T>[g^] exp(i S EH [g^}), (1.3) 

JGeom(M) 

where g'^, g'/j are induced boundary geometries, is called the propagator 
by analogy with Feynman path integral for a single particle and can be re- 
garded as a transition amplitude between an initial and final spatial geometry. 
(Which is illustrated in Fig. 11.11 ) 
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Figure 1.1: A quantum transition between an initial and final spatial geom- 
etry via quantum superposition of intermediate geometries described by the 
gravitational path integral ( 11. 3p 

Unfortunately, there is a number of mathematical problems associated 
with (|1.2|) . (I1.3p . let alone the physical interpretation (cf. [T], j2]). First and 
foremost, as always in quantum field theory one has to take care of conver- 
gence of a functional integral. This can not be circumvented here by applying 
a perturbative expansion, therefore a suitable (non- violating diffeomorphism 
invariance) regularisation and renormalisation of the path integral itself is 
necessary. Secondly, since there is no absolute parametrisation of the space 
of geometries Geom(M) we do not know what we are integrating over. Espe- 
cially, how the diffeomorphism-covariant (because of the integral having to 
be diffeomorphism-invariant) measure P[<7^] should be defined. As a result, 
one has to start with gauge-covariant (i.e. diffeomorphism-covariant in this 
case) fields and fix the gauge afterwards but this gives rise to Faddeev- Popov 
determinants, which are exceptionally difficult to be computed nonpertur- 
batively. Another issue is the very existence of the measure. Thirdly, for 
any choice of variables the action ( II. ip is quadratic in fields, which hinders 
evaluation of integrals in (II. 2p . ( II .3p . Furthermore, integrands of those are 
complex and therefore, a kind of Wick rotation i.e. an analytic continuation 
such that t H- i t is required. Due to the dependence of all metric components 
on time as well as non gauge-invariant nature of the above straightforward 
Wick rotation it is not so obvious what it should look, however. Were all 
these problems tackled, for d > 3 there is yet another, known as the con- 
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formal factor problem. Briefly speaking, it consists in unboundness (from 
below) of the action (II. ip . triggered by the conformal mode of metric. This 
makes the Wick-rotated path integral still potentially divergent. Going be- 
yond this section, in the context of causal dynamical triangulations the issue 
was addressed and probably resolved [1] (see also [2]). 

Last but not least, in the partition function (II. 2p (as well as in (II. 3p ) 
the topology of spacetime M is fixed. In principle, one might consider an 
extension of (11.2)1 including a sum over topologies 



Although this makes the situation even more complicated, it is argued by 
some (see for instance [3]) to be necessary from a point of view of the theory 
of quantum gravity. Nevertheless, there is still no experimental/observational 
evidence for the possibility of topology fluctuations which would arise then. 
In any case, this issue will be left out here. 

1.2 Method of dynamical triangulations 

There are several approaches [I] to fix the first of the above mentioned prob- 
lems with the gravitational path integral i.e. finding a suitable regular isat ion 
and renormalisation. (Which - hopefully - should also help to resolve the 
other issues.) Inspiration has been coming especially from lattice gauge the- 
ory and its successes in quantum chromodynamics. In this way so-called 
quantum Regge calculus appeared [6], originating from Regge calculus [5] 
a method of approximating solutions of classical gravitational theory with 
piecewise linear manifolds. It is characterised by the fact that one introduces 
a lattice regularisation in the form of a discretisation of geometries and the 
path integral is over discretised geometries only. It ought to be noted that 
in contrast with (classical) Regge calculus it is not an arbitrary continuous 
manifold that is to be approximated here but the integral over all such mani- 
folds. Hence, discretisation of geometries does not need to be arbitrarily fine. 
For example, a quite limited number of lattice plaquettes can be considered. 
Nevertheless, it should be taken to infinity along with a plaquette size taken 
to zero in the limit of a presumable renormalised - i.e. continuum - theory. 

Dynamical triangulations (DT) constitute a distinct variant of such a 
scheme (see jlj). They are characterised by the integration over piecewise 
linear manifolds assembled from polytopes (a generalisation of polyhedra) 
which have fixed lengths of the edges, functioning as a gauge- invariant ul- 
traviolet cutoff. Polytopes are traditionally regarded as being intrisically 




Geom(M) 





6 



Minkowskian (i.e. in a sense flat - which makes manifolds piecewise lin- 
ear), although this assumption may require a change for causal dynamical 
triangulations. (Generally, there is no extrinsic geometry, since manifolds 
assembled from polytopes are not embedded anywhere - the same as in the 
case of continuous ones.) 

More explicitly, in the method of dynamical triangulations the continuous 
partition function (11.21) is replaced with the discrete 




where k, A denote coupling constants of the regularised theory, correspond- 
ing to the inverse of gravitational constant G and cosmological constant A 
(divided by G, to be precise), respectively (cf. (11.9P here below); the sum, 
to which the integral reduces, goes over all "polytopialisations" (e.g. trian- 
gulations, cubulations) T of a fixed topology, weighted by the exponential 
of the regularised action S^fT] and normalised by the factor Cj (rank of 
the automorphism group of T). The latter replaces the measure Vlg^} and 
does not allow for overcounting of geometries, which would violate diffeo- 
morphism invariance. Such a description of geometry is stripped of the use 
of coordinates, though in principle a coordinate system might be assigned to 
each polytope. This is absolutely unnecessary, however, because all geomet- 
ric properties are contained in the connectivity of a polytopialisation (the set 
of neighbouring relations between polytopes), as it will be discussed below. 
Furthermore, thus we obtain a gauge-invariant parametrisation of the regu- 
larised Geom(M), a desire for which was expressed in the previous section. 




Figure 1.2: Curvature at a vertex of a 2-dimensional simplicial manifold: 
positive on the left (5 > 0) and negative on the right (5 < 0); the general 
case, with non-equilateral triangles 
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Since ( II. 5p should be (and - as it was e.g. explicitly shown [7] for d = 2 - 
is indeed) independent of a polytopialisation usually the simplest and most 
convient choice is made by taking only simplicial manifolds (i.e. built from 
simplices), called triangulations, into consideration. Moreover, the simplices 
are assumed to be equilateral (in Minkowskian sense). This is exactly what 
will be done in all the following. Hence the symbol T above as well as 
Sr, denoting the Regge action (see (jl.9p here below), and the very name of 
dynamical triangulations. 

The detailed way of encoding of geometry in a simplicial manifold is as 
follows (cf. [2], [8]). Actually, it is enough to consider curvature alone. In 
general, for d = 2 sectional curvature on a manifold coincides with so-called 
Gaussian curvature and Ricci scalar is equal 2 times the latter. The Gaussian 
curvature can be determined from the parallel transport of a vector around 
a closed loop. On a simplicial manifold such a loop has to circle a vertex, 
otherwise no effect will be ever observed (since simplices are flat). A vector 
will then change its orientation if and only if a sum of the adjacent angles 
^ 27T. A difference is called the deficit angle S. This is exemplified in Fig ll.2l 
For d > 2 one has to sum over all 2-dimensional submanifolds dual to a certain 
(d— 2)-dimensional subsimplex to obtain the nontrivial Ricci scalar, localised 
at this subsimplex. Hence, curvature is associated with (d — 2)-subsimplices 
and expressed by deficit angles 5 around them, which - by generalising a 
notion for d = 2 - are defined through sums 



where 9^ are ditopial (a generalisation of dihedral) angles adjacent to a given 
{d — 2)-subsimplex A^ -2 ^. Deficit angles, as can be observed, are determined 
by the manner in which simplices are glued together - their connectivity. 

From the above discussion, there is the following correspondence between 
terms of the non- regularised (formula fll.lH ) and regularised Einstein-Hilbert 
action 



where vol denotes the volume of a simplex/subsimplex (cf. Section 1.4). 
Consequently, the Regge action amounts to 
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Its precise form for d — 4 will be given in Section 1.4 (in the Wick-rotated 
version, yet for the modified regularisation, as will become clear soon). 

1.3 EDT and their failure 

The model whose regularisation scheme was outlined in the previous section 
is called Euclidean dynamical triangulations (EDT). Unfortunately, there are 
severe problems (cf. [2], [8]) with moving further in EDT, towards a certain 
nontrivial and physically relevant continuum limit, which is the ultimate aim 
of the framework. Everything starts with the absence of a prescription for 
Wick rotation of the partition function (II. 5p . 

where Sr\T'\ is the Euclidean (Wick-rotated) action. Therefore, the ad hoc 
Wick rotation is applied, in which Lorentzian geometries are replaced with 
Euclidean ones just by hand. As a result, the summation has to go over 
simplicial manifolds devoid of the notion of time and causal structure. It 
is not known how to retrieve the latter afterwards. Almost certainly, it can 
not be done by straightforward inverse of the Wick rotation, also because 
of the occurrence of Euclidean geometries anomalous from the perspective 
of causality. Namely, the problem consists in changing of the topology [5] 
between (d — l)-dimensional submanifolds, especially by events of splitting 
of such a submanifold into at least two disconnected ones, as it is exempli- 
fied in Fig. 11.31 For a Lorentzian manifold this amounts to disintegrating of 
space during its time evolution, in other words, branching of spacetime into 
so-called baby universes. (Obviously, the time inverse of the process in the 
form of joining of the disconnected space components should be identically 
likely yet it is excluded by the requirement of the fixed spacetime topology.) 
A consequence is the degeneration of light cones; causality is violated because 
at a branching point there is more than one possibility for further course of a 
test particle's world line. Moreover, the above phenomenon turns out to be 
generic and highly degenerate geometries are dominating in the Wick-rotated 
(formula (11.101) ) counterpart of the sum (11. 5ft . 

Hence, we are stuck in the Euclidean theory and this is where the very 
name of EDT derives from. Bad news has been confirmed by results (see 
[1], [2]) from numerical simulations. Namely, there is no relevant continuum 
limit for d > 3. Rather than that, the space of coupling constants k, A for 
the model of d = 3, 4 (and probably any d > 3) contains a critical line X cr it{k) 
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(where the continuum limit exists), at which there are two extreme, unphys- 
ical phases, separated by a bicritical point. A useful tool to study them is a 
standard notion of (effective - in the continuum limit sense) dimensionality, 
called Hausdorff dimension, which works well also on discretised manifolds. 
In the phase of small fc's, called the crumpled phase, simplicial spacetimes 
(or rather Euclidean spaces) of a tiny size are produced. Moreover, there 
exist several vertices to which almost all simplices are connected. Therefore 
distances between the latter are small and in the continuum Hausdorff di- 
mension of space is going to oo. On the contrary, in the phase of large k's, 
known as the branched polymer phase, generation of baby universes is the 
dominating process. Therefore every simplex has only a limited number of 
neighbours and in the continuum Hausdorff dimension approaches 2, which 
enhances interpretation of spaces in this phase as a particular kind of frac- 
tals - namely branched polymers. Neither of the phases yields Hausdorff 
dimension equal or at least close to the expected d. 




Figure 1.3: A change of spatial topology by branching of spacetime with the 
resulting behaviour of light cones shown 



1.4 Imposing causality on DT 

Causal (or Lorentzian) dynamical triangulations (CDT) were invented [H] 
as a response to the failure of EDT and constitute their modification [10] , 
[TT] . The essential idea behind this approach is to consider from the outset 
only causally well-behaving simplicial manifolds. More precisely, a triangu- 
lation must admit a global proper-time foliation. It is only when it is built 
from layers of simplices to which a proper time variable can be assigned - 
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hence space globally evolves in time (spacetime is globally hiperbolic). Con- 
sequently, the topology of space is conserved, which means the exclusion of 
spacetimes branching into baby universes. However, it ought to be stressed 
that changes of spatial topology are not excluded from consideration in prin- 
ciple. Individual geometries in the integral fll.2p /sum (ll.5p are in general 
highly nonclassical, as it is always in quantum physics, and possibly could 
also include fluctuations of spatial topology. Causality which would be lost 
in such a case might as well be restored for a superposition of quantum space- 
times, corresponding to a classical solution of gravitational theory. On the 
other hand, we do not even know whether the restriction on topology sur- 
vives in the continuum limit of the CDT model. Nevertheless, the failure of 
EDT has shown that in the method of dynamical triangulations variability 
of spatial topology has to be treated exceptionally (see [8]) if at all (and its 
total rejection is worth trying). The same concerns variability of topology of 
the entire spacetime, discussed at the end of Section 1.1, which leads to the 
super-exponential growth of the number of geometries in (11.51) . 





Figure 1.4: Types of simplicial building blocks in CDT (without the dual 
ones) for d = 3 (up) and d — 4 (bottom) 

Due to the layered structure of simplicial manifolds in CDT there is a 
limited number of allowed building blocks from which they are constructed - 
simplices placed differently with respect to the time direction. Specifying the 
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statement from the previous paragraph, a proper time variable is assigned 
to simplicial vertices and (d— l)-dimensional surfaces spanned between such 
at a given time (which obviously are forbidden to intersect). As can be seen, 
time is discrete. As a result, we have d types of simplicial building blocks, 
denoted by numbers of their vertices lying at times t and t + 1 (where t is 
arbitrary). In the most interesting, lowest dimensions we have the following 
ones. For d = 2 there are simplices of types (2, 1) (2 vertices at a time t, 
joined by a spatial link and one at t + 1; the other two links are timelike) 
and (1, 2), which is just time inverse of (2, 1) (I call them dual to each other). 
Types for d — 3, 4 are depicted in Fig. 11.41 For d = 3 there is a pair of dual 
types of simplices, (3, 1), (1, 3) (with a spatial triangle) and type (2, 2) (with 
two spatial links only). For d — 4 there are two dual pairs: (4, 1), (1, 4) (with 
a spatial tetrahedron) and (3, 2), (2, 3) (with a spatial triangle and a spatial 
link). In all the following I will be denoting collectively the latter as {4,1) 
and (3,2), respectively. 

The next step in construction of the CDT model (discussed extensively 
in [TS]) is the introduction of anisotropy between time and space. Squared 
lengths of timelike and spacelike simplicial links are set as, respectively, 

7 2 — — 2 7 2 — 2 /i ii\ 

'timelike — Oi a , { space iik e — d , \ 1 • 1 1 ) 

where a is the cutoff size and a > the anisotropy parameter. Hence sim- 
plices are generally no longer equilateral in the sense of the Wick-rotated 
model (in EDT we had a = 1). Now according to the correspondence (11.71) . 
( 11.81) one has to calculate volumes of timelike and spacelike simplices and sub- 
simplices as well as values of ditopial angles - all as functions of a. (Different 
conventions may be used.) Irrespective of what might seem, all Lorentzian 
volumes are real and > but angles in general are complex. Yet it is not a 
problem, as we eventually obtain a real regularised action. Its form without 
substituted volumes and angles (for simplicity) - an incarnation of (II .9p - 
for d — 4 is 

S R [T] = k £ vol(Ai 5 <°)) 5 -f + k vol(A^H- 

A £ vol(A^))-A ]T vol(A^), (1.12) 

A u>i) cT Af\r 

where by (3, 0) I denoted spacelike triangles and by (2, 1) timelike ones. 
What is essential is that the action (I1.12p can be naturally Wick-rotated. 
The analytic continuation, such that a y —a, goes through the lower half 
of the complex plane and in the lowest dimensions can be defined: for d = 2 
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always, for d — 3 if a > |, for d — 4 if a > X. The bounds on a originate 
from the form of ditopial angles. Next, we perform sums in (11.121) and use 
topological identities (by Dehn and Sommerville) , such as the Euler identity, 
for the numbers of simplices and subsimplices of all types. Finally, for d = 4 
and the assumed spacetime topology S 1 x E (the periodic boundary condition 
in proper time) where £ is 3-dimensional and compact we get the Regge 
action which has the Euclidean form 

Sr(T) = -(ko + 6A)iV + k,(n¥' 1] + Ni 3 ' 2} ) + 

A(2N { 4 4 ' 1] + iVf i2) ), (1.13) 

( A I s ) (3 2} 

where Nq denotes the number of O-simplices (vertices), N\ , AQ the 
numbers of 4-simplices of the corresponding type and Ko, k 4 , A are coupling 
constants and quite complicated functions of k, A and a. However, their 
dominating behaviour is determined by single variables - k, A and a, respec- 
tively, and in particular A = 4^ a — 1. Consequently, the overall form of 
the (Wick-rotated here) partition function (II. 5p is slightly altered into 

Z(k ,k,,A) := V 7 ^exp(-5|[T]), C r = |Aut(T)|. (1.14) 
r ° r 

Phase structure of the CDT model for d = 3, 4, which can be inferred from 
numerical simulations (whose technique will be brought up in Chapter 2), 
differs widely from that of EDT. Actually, there is a disparate continuum 
limit already for d = 2 [5] but I will not discuss this fact here. Three coupling 
constants instead of two amount to the 3-dimensional space, containing a 
critical surface ^^(kq, A) (where the continuum can be approached) i.e. the 
phase space of the model. As Fig. 11.51 indicates, we observe [13], [H] three 
phases: A, B and C. They have the following characteristics. Spacetimes 
in phase A (of large re's) consist of numerous small, disconnected universes, 
lying along the time direction and as such they are physically irrelevant. 
This phase is considered a remnant of the branched polymer phase from 
EDT. Spacetimes in phase B (of small A's) are no more useful, being very 
short in time and having at each temporal end a vertex of a very large order 
- the same as in the crumpled phase from EDT. Therefore phase B can be 
regarded as a remnant of the latter. Fortunately, there is an additional phase 
C, in which relevant spacetimes - extending both in time and space - are 
observed. Its existence has been a promising result of the CDT model. It was 
shown [15] that a universe from phase C has indeed Hausdorff dimension d. 
Furthermore, recently it turned out [T4] the above phase structure is related 
to so-called Horava-Lifshitz gravity. 
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Figure 1.5: Phase space of CDT for d — 4 with the measured boundaries 
between phases A, B, C (from [H]) 

One remark has to be added. All along, from the very beginning we are 
dealing with gravitational field in vacuum i.e. without any matter fields. 
This is the work on pure quantum gravity (or quantum geometry). However, 
the inclusion of other fields in the gravitational path integral has always been 
considered and also in CDT is subject to ongoing research. 
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Chapter 2 



Monte Carlo simulations and 
the semiclassical solution 

2.1 Numerical setup of CDT 

To solve the model of causal dynamical triangulations, as it is formulated 
in the previous chapter, one has to compute the sum (jl.l4)l . Thus far, it 
was carried out analytically only for d — 2 |9j (see also [11]) and there is 
little hope for moving towards higher dimensions. However, like every lattice 
gauge theory CDT (as well as EDT) is well-suited to be investigated numer- 
ically, provided one can ensure the regularised version of space Geom(M) is 
sampled thoroughly in simulations. Hence in this context the problem with 
the regularisation may be described as twofold: not only has the space of all 
possible triangulations to be appropriately dense in the Geom(M) but also 
numerical simulations enable us to approach every element (in the certain 
class) of the latter. The first issue is so far assumed, the same as the fact 
that £p in (OP is the measure. To address the second one, Monte Carlo 
methods are applied [12]. More specifically, to approximate fll,14p we gener- 
ate some number of simplicial manifolds, which are obtained consecutively 
from the previous ones - starting from a certain ad hoc constructed, very 
simple triangulation - by local transformations called moves. Factor is 
automatically included there. Now, the moves are hoped to be ergodic in the 
sense of the possibility of reaching every triangulation from any other - in a 
certain class - by the successive application of them. What I mean by the 
"certain class" will become clear below. Concerning Monte Carlo moves, they 
consist in replacing of a complex built from some number of simplices with 
another complex, in general constructed from a different number of building 
blocks yet respecting the topological and causal structure of a triangulation, 
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in conformity with the assumptions of CDT. For d = 2 there is a set of 3 
such moves, for d = 3 there are 5 of them and for d = 4 there are 7 - their 
detailed description (omitting d = 2) is contained in [12J . 

From now on the discussion will be restricted to d = 4, which I actually 
studied. Regarding the argument from Section 1.4, the topology in the reg- 
ularised path integral (11. 14j) is fixed. For convenience one may consider a 
topology with a periodic boundary condition in proper time, S 1 x £ where E 
is 3-dimensional and compact, as it is in the case of the action f II . 1 3 j) . This 
is exactly what we do in the discussed Monte Carlo simulations, choosing 
£ = S 3 (thus far - as the simplest case). Hence the topology is S 1 x S 3 , 
which is the first constraint characterising the class of triangulations which 
I investigated, refering to what was said in the previous paragraph. Another 
possible condition is the fixed number of simplicial temporal layers, in other 
words, the number of discrete proper times i, which in the discussed Monte 
Carlo simulations is set as T = 80. The parametrisation of time I will be 
using here is such that % G [0, 79] and i = T = 80 = 2 = (the periodic 
boundary condition). Furthermore, with the aim of keeping the number of 
4-simplices in a triangulation roughly within a certain range, in the case 
considered here we modify the action (11.131) . 

(T) S$(T) +e((Ni^) -M) 2 , (2.1) 

where (AQ 4 ' ^) is the average total number of 4-simplices of type (4, 1) (see 
the next section), e is a small constant and A/4 = 80000 (some other values 
were also studied - as one can check in [15], [16], [T7], [18] - yet for the 
modification of the action different from (12. ip ). It ought to be noted that in 
the considered topology the ivj 1 ' 4 ^ is exactly the same as the A^ 4 ' 1} , because 
only those types of 4-simplices contain spatial tetrahedra (cf. Fig. II. 4p and 
hence every (4, 1) simplex has its counterpart (1,4), sharing a spatial tetra- 

(3 2) 

hedron with it. On the other hand, the AQ is by construction related to 
the ivj 1 ' 4 ^ . Thus the number of all 4-simplices fluctuates in a certain range of 
values. Last but not least, concerning the discussion of the phase structure 
of CDT presented in Section 1.4, in order to find ourselves at the critical 
surface we have to fine-tune k^. This is done here using the requirement 
(ivj 4 ' 1 ^) = A4 i.e. the average number of (4, 1) (as well as (1,4)) simplices is 
kept fixed. Consequently, e characterises fluctuations of (A 7 "] 4 ' 1 ^) around A/4. 

All what follows concerns numerical results from the Monte Carlo simu- 
lations carried out in the physically relevant phase C. I conducted the most 
extensive analysis for values of the coupling constants k = 2.2, A = 0.6 but 
all what will be said in this and subsequent chapters is true at least for a cer- 
tain range of k,q and A in phase C. All of the presented graphs will be given 
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for e = 0.00004 yet e = 0.00008, 0.00002 alter them only slightly. Further- 
more, the results are practically the same also for the variant of ( 12. ip with 
(A/" 4 ) (i.e. all simplices) instead of (AQ ), with s = 0.00008. The number of 
statistically independent (which means separated by a sufficient number of 
moves) Monte Carlo configurations (i.e. triangulations) generated in each of 
the above cases was about 10 5 . 

2.2 The semiclassical limit 

From now on I will be using the notation N p>q = A/4 (and blank indices will 
mean an undefined type of simplices). The most fundamental observable in 
numerical simulations of CDT is the average number of simplices of a given 
type at a given time 



where n is the number of configurations we average over. For the already 
mentioned total number of (p, q) simplices (N Ptg ) one simply has to sum 
over all times t. The above expression is an estimate of a grand canonical 
observable for the system fluctuating around (A^i) = A/4, provided all con- 
figurations are typical i.e. of large probability. Actually, calculating (12. 2p 
one has to take one more thing into account. (I restrict myself to (4, 1) 
simplices at the moment - which is the same as restricting to (1,4) - but 
cf. Section 3.1 and Appendix A.) Namely, looking at an arbitrary temporal 
distribution N^(i) in phase C we observe two regions: one where a universe 
is well extended in spatial directions and one where it is almost not at all 
(since (4, 1) simplices correspond to spatial tetrahedra). In other words, the 
extended universe is localised in a certain range of discrete times i. Since 
CDT is a priori symmetric with respect to translation in the time direction, 
we reparametrise time for every triangulation in such a way that the so-called 
center of volume of N^i(i) lies at % — 39. (However, as it will be seen below, 
the effective center of volume is between i = 39 and i = 40, at t = 39 + |.) 
More precisely, the center of volume icy is defined here as %' for which 



reaches a minimum; if there is more than one, A^i^cv) must be the greatest 
among them. (There are alternatives to (12. 3p but they all yield the same re- 
sults up to a difference of 1.) Now we may calculate the average distribution 




(2.2) 



k=l 



(* + I) #4,1 (1 + m ° d ( T + *" + i - 1, T)) 

1 2 



(2.3) 
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(N^i(i)). The result is shown in Fig. 12.14 along with the average fluctuations 
of N^i(i) (see Section 3.1). One can distinguish three regions in time evolu- 
tion: the "stalk", where the spatial extension of the universe is very small, 
consisting of only several spatial tetrahedra (but not less than 5, which is the 
minimal number for a simplicial S 3 ), the "blob", where the universe is well 
spatially-extended (with the double maximum at % — 39, 40) and the interme- 
diate "tail" (see below). Moreover, the calculated (iV 41 ) = A/4 up to a very 
small deviation. It has yet to be added that (N^i(i)) should be invariant 
under time inversion. Therefore, it is additionally symmetrised with respect 
to the effective maximum t — 39 + ~, which corresponds just to doubling of 
the number of configurations we average over (and this is justified by the fact 
(-^4,1 (0) i s very symmetric even without that). 




Figure 2.1: The average distribution (N^^i)) with the semiclassical solution 
fitted and the average quantum fluctuations (see Section 3.1) denoted by 
vertical lines (not to be confused with errors, which are invisible here!) 

It turns out [T7] that the blob region of our spacetime can be fitted sur- 
prisingly well with a function 

N cos 3 {u{t-t )), (2.4) 

where to = 39 + | and No, u are certain constants. This is illustrated in 
Fig. 12.11 As one can see, at least in magnification, the fit starts to deteriorate 
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close to the stalk and this is one of the symptoms of being in the transition 
region called the tail above. Furthermore, taking the dependence on A/4 into 
consideration, (I2.4p has the form with explicit scaling 

A s JVJ cos 3 {A t AfP(t- t )), (2.5) 

where A t , A s are certain constants. Hence it is a limit of the CDT model, 
an effective solution. What is crucial in this statement is that eventually 
the limit may be called semiclassical [TH] - the semiclassical solution of CDT 
- because it is [17], [32] the (classical) solution of the Wick- rotated mini- 
superspace model, whose action 

Sco„t = ^jdt v^(^|^ + ^(tfs - AV 8 (f)) , (2-6) 

where V3 (t) = 27r 2 a(t) 3 (where a(t) denotes the scale factor of a homogeneous, 
isotropic universe) is the spatial 3-volume, ki = 9(27r 2 )3 and A is a Lagrange 
multiplier fixing 

Jdt ^~ t V 3 (t) = V 4 , (2.7) 

where V4 is a certain 4-volume of spacetime. The solution of ( 12 .6p is a ge- 
ometrical 4-sphere, which amounts to the Euclidean de Sitter universe, the 
maximally symmetric spacetime in the case of a positive cosmological con- 
stant. This comes as a great surprise, since individual triangulations - as it 
was stressed before - are highly nonclassical and a priori devoid of any sym- 
metry, except the one with respect to a time translation. What is even more 
important is that also dynamics of the model i.e. in particular (quantum) 
fluctuations of iV 4il (i) can be described in that framework - by the effective, 
semiclassical action (12. 6p itself [TH]. This is what I will deal with - to some 
extent - in the remaining chapters. 
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Chapter 3 



The single slice structure 



3.1 Expansion of the semiclassical action 

By construction, simplicial building blocks in temporal layers of a triangula- 
tion have to be connected in an appropriate order. Namely, in the simplest 
case, between a (4, 1) simplex at a time % and another one at % + 1 one has 
to put consecutively a (3,2), a (2,3) and a (1,4), with successive simplices 
sharing a face (i.e. a tetrahedron), as it was already said for the case of (1, 4) 
and (4, 1) simplices, which share a spatial tetrahedron. (Generally, each of 
the above mentioned five simplices may occur multiple times there.) Hence 
the layer has some internal structure. However, in this chapter I will restrict 
the discussion to spatial hypersurfaces formed by common faces of (4, 1) and 
(1,4) building blocks - forgetting for a while about the other types of sim- 
plices - in accordance with Section 2.2. This is what I call the single slice 
structure: the structure in which a temporal layer is contained between two 
spatial slices (hypersurfaces) at times i and i + 1, the second of which be- 
longs already to the next layer and there is nothing else, (3, 2) (alternatively 

(3.2) and (2,3)) simplices are in some sense integrated out. As it will be 
shown in the next chapter and sketched out in Appendix A, actually there 
is a nontrivial structure emerging from (3, 2) simplices, as well as (3, 2) and 

(2.3) considered separately, compatible with that of (4,1). 

According to Section 2.2, the average distribution (A^ 41 (z)) is well de- 
scribed in the blob by the semiclassical solution (12.51) . resulting from the 
classical, continuous action (12. 6p . Since time i is discrete (N^iii)) actually 
corresponds to a discrete version of (12.51) . Therefore we expect that in the 
analysed semiclassical limit CDT can be described (at least in the blob) by 
the discretised mini-superspace action S^ s , whose explicit form will be dis- 
cussed in the next section. More specifically, to take the quantum effects - 
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fluctuations around the effective, semiclassical geometry - into account we 
make the expansion of around the distribution N^i(t) = (iV^^i)), which 
should extremise the action, with an increment SN^i(t): 

S dis {N 4A (t) + SN^t)] = S dis [N 4A (t)}+ 

1 T_1 

2! S 9 NiM3) d N4Al) S dis [N AA (t)} 5Nt A (i)5N AA (j) + 0(5iV 4il (t) 3 ), (3.1) 

' j,i=0 

which in particular implies 

Vz: d N4tl[i) S Ais [N^(t)] = 0, (3.2) 

vanishing of first derivative of the action. 

The second most fundamental observable in our Monte Carlo simulations 
of CDT is the covariance matrix of the number of simplices of a given type, 
whose element - the correlation between given times t' and t" - is defined as 

C P , q (t',t") = i£(tf«(0 - (tf«(0» (JV$(0 - {Ng>(?))) (3.3) 
fc=i 



(cf. 12.21) . Quantities ^C P)q {t, t), included in Fig. 12.11 for (p,q) = (4,1), are 
the average fluctuations around the expectation values (N p ^ q (t)). Now, let us 
again restrict ourselves to (4,1) simplices. Analogically with (AT 4)1 (t)), 
are centered at icv' s an d is additionally - but on the justified basis - 
symmetrised with respect to t = 39 + |, which corresponds to doubling of the 
number of configurations. In our numerical setup matrix is invertible 
and hence the inverse (of the) covariance matrix 

P 4>1 = Cll (3.4) 

can be found. What is essential is that there is the relation 

(Pi,l)i,j = ^iv 4jl ( : ,-)9Ar 4il (i)5'dis[iV4 j i(t)] ) (3.5) 

where I introduce the notation i,j = i + l,j + l for matrix elements (since 
time G [0, T — 1]). As a result, from analysis of matrix P 41 one can infere, 
at least to some extent, the exact form of the action S^is- It turns out the 
inverse covariance matrix has only simple, tridiagonal structure 



P 



4,1 



( a,\ bi * * br \ 

bi CL2 ' ■ * * 

* ' • . ' • . ' • . * 

* * ■ • olt—\ bx— i 
V bx * * bx-i clt ) 



(3.6) 
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with two additional elements at the antidiagonal resulting from the periodic 
boundary condition in time (neglecting numerical noise, denoted by the as- 
terisks) and Vi: aj > 0, Vi: 6j < 0. It is also by construction symmetric with 
respect to both the diagonal and antidiagonal (a^ = ax+i-i, h = &T-mod(i,T))- 



3.2 Form of the discrete action 

The discretised action, with the modification (12. ip taken into account, has 
the general form 

T-l 

S dis = h ]T (S k (t) + S p (t)) + e «iV 4) i) - M) 2 , (3.7) 

i=0 

where kij^t^kit) = &i<Sfc[-W4,i(£)], ^i E< ^pW = kiS p [N4 t i(t)] are respec- 
tively its kinetic and potential part, corresponding to the terms 9 y^ and 

k 2 V 3 (t)3 — \ V 3 (t) in the continuous action (12. 6 p and k x denotes a coupling con- 
stant proportional to the inverse of G. In all the following I call Sk[N^i(t)], 
SplN^iit)} the kinetic and potential term, respectively. I assume that the 
potential term S p (t) = S p (N^i(t)) and has the form 

S p (t) = hN^irf - XN^t), V(t) = hS p (t), (3.8) 

where k 2 , A are coupling constants proportional to k 2 and A, respectively. The 
kinetic term is more involved. However, from (13. 5p . (13. 6p one immediately 
infers that §k(t) = Sk(N^i(t), A^ 41 (t + 1)). Consequently, the diagonal of P 4j i 
depends on both the kinetic and potential part of the action Sdi S , whereas its 
super-/subdiagonal on the kinetic one alone. I assume that the kinetic term 

s (f) = (N iA (t + I) - N 4A (t))* (jy 4|1 (t + l)-^ 4|1 (t)) 2 
feU h(N 4>1 (t),N 4A (t+l)) h(N 4A (t+l),N 4jl (t)) { - ] 

and satisfies the relation 

dN^i-vd^SklN^it)) = 0, (3.10) 

which obviously happens at least for every h(N 4> i(t), N 4) i(t + 1)) that is a 
quotient of combinations of power functions. According to (13. 5 p I may define 
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the analogue of (13.101) : 

#4,1(1) #4,1 W 

which - keeping in mind the form of (13.81) and (13. 9p - in the relevant re- 
gion i.e. at least in the blob, should be equivalent to the non-eliminated 
kid^s^SplN^i^b)] = — lkik 2 N^i{t)~3 . Hence I can check the potential 

term (13. 8p independently of the explicit form of Sfc(i) i.e. form of the func- 
tion h(N^x(t), N4 t i(t + 1)) (Fig. I3.ip . (In figures I use the abbreviated no- 
tation di = 9jv 41 (i), [t] = [N^i (t)).) The obtained V{i) exhibits in the blob 
the considerable numerical noise. As it will become more evident below, the 
underlying reason is a small magnitude of the potential term. 



O 



0.0001 



le-05 



le-06 




log[<AMi)}] 

Figure 3.1: — |fcA^ 4)1 (i)"i fitted to V(z), in the relevant region where V(i) < 



The simplest form of the kinetic term (13. 9p one can consider is 
~ 2(iV 4 , 1 (t + l)-iV M (t)) 2 

Skit) = N^ + N^t + l) ■ (3 ' 12) 

Now, regarding the form of the semiclassical action (13. 7p . the extremum 
condition (13.21) in the expansion (13. ip yields 

Vi: -^ 4il(i) ^[JV 4 ,i(*)] ^d N ^S p [N Ail (t)] + ^e((N 4tl )-M A ) (3.13) 
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Figure 3.2: \k 2 N^i{i)'~^ ~ ^ fitted to — dj^^^SklN^^t)]; for the illustrative 
purposes the term with e is absorbed into A 



and I may examine the validity of (13.1 1) . (13.21) by fitting first derivative of 
the potential term (13. 8p to first derivative of the kinetic one (I3.12p (Fig. 13.21) . 
This clearly does not work in the stalk. In the blob and perhaps the tail, 
making the use of (13.51) . for the superdiagonal of it should be 

Vz: h d N4 l{i+1) d N4A{i) S k [N 4A (t)} = (P 4 ,i)m+i - 2e (3.14) 

(a choice of the subdiagonal would obviously imply the same). However, 
the kinetic term (I3.12p starts to fail there. I analysed different functions 
h(N 4tl (t), N 4A (t + 1)) and something like \ (y/N Ail (i) + y/N 4>1 (i + l)) 2 turns 
out to behave much better. Alternatively, one may consider a modification 
of (I3.12p in the form of a series with coefficients 

~ (iV M (t+l)-iV M (t)) 2 / 

fl k(t) + ^ + i)) (3 - 15) 

Thus far, our numerical precision is sufficient to find only £i and hence con- 
straints Vj > 1 : £j = are assumed. (£1 = 2 corresponds then to an approx- 
imation of \(y/N^i(i) + \/N^i(i + l)) 2 and £1 = is obviously equivalent to 



24 




100 1000 



log[i((AMi)) + <AMz + l)))] 
Figure 3.3: fcj d N4 ^ i+1 )d N41 ^S k [N Ajl (t)) fitted to (P 4 ,i)»,»+i - 2e 
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log[(^4,iW}] 

Figure 3.4: -pJV^^-f fitted to (P 4 ,i)m - 2e - A* ^ l(i) 5 fc [iV4,i(t)], in 
relevant region where the latter < 
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(|3.12p .) I have examined (I3.14p with the above kinetic term fitting simultane- 
ously ki and £1 (Fig. 13 .3j) and obtained a value of the latter ^ = 0.327±0.030, 
closer to \ than to 0. Concerning the relation (I3.13p . presented in Fig. 13. 2\ 
its behaviour is practically unchanged by the modified kinetic term (I3.15p . 
Therefore, I do not show the corresponding version of the graph. At last, 
taking the values of k\ and £i from Fig. 13.31 1 can investigate the potential 
term (13. 8p at the diagonal of P 4j i, since 

Vz: hdl iA{) S p [N 4>1 (t)} = (P,, 1 ) i , i -2e-k 1 d% 4 ^ ) S k [N 4>1 (t)}, (3.16) 

which results from (13. 5p (Fig. I3.4p . Now, I may compare the values of the 
constant k = k\k 2 obtained in Fig. 13.11 and Fig. 13.41 They amount respec- 
tively to k = 0.908 ± 0.033 and k = 0.873 ± 0.021. Furthermore, I may 
compare the value of k 2 from Fig. 13.21 (or rather its variant with the modified 
kinetic term) with the one that I can calculate from Fig. 13.41 They amount 
respectively to k 2 = 18.9 ± 1.1 and k 2 = 22.97 ± 0.81. The presented uncer- 
tainties include only the accuracy of fits and do not contain estimates of the 
systematic errors. 
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Chapter 4 

The double slice structure 



4.1 Evidence for the double structure 

According to Section 3.1, (3,2) simplices bond together the adjacent tem- 
poral layers of (4,1) building blocks. Furthermore, refering to Fig. I1.4[ one 
can observe that (4,1) simplices contain - apart from a spatial tetrahedron 
- only (3, 1) timelike tetrahedra, while (3, 2) simplices possess both (3, 1) 
(two) and (2,2) (three) tetrahedra. Hence every (3,2) simplex has to share 
a face with three others. Numerical simulations confirm that this type of 
building blocks forms closed hypersurfaces in the sense of the dual lattice 
and thus we may speak about layers of (3,2) simplices. Therefore I extend 
my description to what I call the double slice structure - characterised by 
the inclusion of remaining (3,2) simplicial building blocks, which implies 
additional spatial slices between those of the (4,1) i-e. an effective second 
slice in every temporal layer, being a spatial section of the (3,2). Neverthe- 
less, either (3, 2) or (2, 3) building blocks are still in some sense integrated 
out, as both types are considered collectively. As it was already mentioned, 
Appendix A makes the outline of their separate treatment. 

I begin with a computation of the average temporal distribution of the 
number of (3,2) simplices, using formula (12.21) . In agreement with what 
one might expect, an arbitrary distribution Ng \(i) turns out to have two 
regions: the spatialy extended and spatialy unextended. Therefore, accord- 
ing to the case of (4, 1) simplices, it requires a translation in time of the 
" center of volume" . However, since I investigate the double structure, distri- 
butions N^l(i) are to be reparametrised also and obviously with respect to 
the same centers of volume. I have found that icv's (defined in Section 2.2) 
for (4, 1) and (3,2) building blocks are identical or at most different by 1 
and hence decided to consider zcv's found for (4, 1) simplices alone. I obtain 
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the average distribution (N 3j2 (i)) having the same three regions as (N^^i)): 
the stalk (with 12 simplices at least, which results from the minimal number 
of 5 (4, 1) building blocks), the blob (with the single maximum at % — 39) 
and the tail. Hence I assume (Ns i2 (i)) = (N 3j 2(i + |)) and fulfilling the 
condition of invariance under time inversion I symmetrise the distribution 
with respect to t = 39 + ~, which corresponds to doubling of the number of 
configurations. Now, it turns out in the blob {N 3j2 (i + |)) is proportional to 
(N^i(i)) in the interpolated sense i.e. I can find a scaling factor p such that 
P (N 3)2 (i + -)) = (N^i(i + !)), where (N^i(i + |)) denotes the interpolating 
distribution i.e. (N^i(i)) symmetrised with respect to i = 39. The joint 
distribution of (iV 4jl (i)) and p (N 3)2 (i + |)) is shown in Fig. 14. 1[ along with 
the average fluctuations, also appropriately rescaled in the case of (3,2) sim- 
plices. Furthermore, the smoothed distribution can automatically be fitted 
with the semiclassical solution (I2.5P and this is done indeed in Fig. I4.1[ up to 
even smaller deviation than in Fig. 12.11 The average total number of (3,2) 
simplices (N 3 ^) ~ 207000. Consequently, the average total number of all 
simplices (both (4,1) and (3, 2)) (N) w 367000 (cf. Section 2.1). 




10 20 30 40 50 60 70 

t 



Figure 4.1: The joint average distribution of {N^i(i)), p (N 3 ^ 2 (i + \)) with the 
semiclassical solution fitted and the average quantum fluctuations denoted 
by vertical lines (not to be confused with errors, which are invisible here!) 
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Using formula (13. 3p one may also calculate the covariance matrix of 
(3, 2) simplices number C$ a and in particular obtain the average fluctu- 
ations around (N 3 ^{i + |)), which are included above in Fig. 14.11 Further, 
I find the inverse covariance matrix Ps^ = C 32 an d it turns out to have 
the same simple structure as (cf. fl3.6j> ) . enhancing the notion of the 
double slice structure. However, the analysis such as in Section 3.2 will not 
be carried out here. Rather than that, in the next section I will focus on the 
joint, double structure of (4, 1) and (3, 2) simplices. 



4.2 The double discrete action 

Since the semiclassical solution (12.51) works well in the blob also in the case 
of the double slice structure I expect that the expansion (13.11) (as well as 
the condition (13.21) ) remains valid for the extended discretised action S^ 1 ^, 
which includes the dynamics of (3,2) simplices. More specifically, 



T-l 



^[ E E E d NrA^N p ^)S^ [N{t)\ SN p , q (i)8N r>s (j) + 0(5N(tf) (4.1) 

' (r,s) (p,q) j,i=0 

and 

^■■d NpAl) S^ ) [N(t)]=0 } (4.2) 

where N(t) = (iV 4 ,i(t), N s , s (t + §)), 5N(t) = (5N 4A (t),6N 3>2 (t + §)) and 
(p,q), (r,s) go over the set {(4, 1), (3, 2)} i.e. in particular there are four 
combinations of derivatives Vi, j in the second term of the expansion. 

The corresponding "double" covariance matrix is defined as the joint co- 
variance matrix of both types of simplices i.e. (13 .3p is extended to 



C dbl (*',0 = -E " <*®0» (N$(f) - (JVg>(0»- (4-3) 



n 

k=l 



Hence Cdbi niay be arranged into four blocks, two of which are the (" single" ) 
covariance matrices 6*4,1, Cs,b an d t wo m i x both types of simplices and which 
will be denoted by superscripts 11, 22, 12 and 21, respectively. (Obviously 
(7<ibi = (Cj^) 1 ".) Each of them is appropriately symmetrised with respect to 
t = 39 + \ . Matrix Cdbi is invertible and the double inverse covariance matrix 
-Pdbi = Cdbi- Now, (13.51) generalises to 

= d NrAj) d NpAl) S^ l) [N(t)], (4.4) 
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where blocks of Pdbh denoted by superscripts, correspond to those of C^u- 
Like P41 (cf. f)3.6p ). the double inverse covariance matrix is very simple 
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(4.5) 



with blocks having band structure (neglecting numerical noise, denoted by 
the asterisks) - except additional elements at their antidiagonals, resulting 
from the periodic boundary condition in time - and Vi : c/ dj Ci > 0, Vz : /, < 0. 
Blocks Plli, Pjbi are also by construction symmetric with respect to both the 
diagonal and antidiagonal and the entire Pdbi is symmetric with respect to 



the diagonal (c/e* = c/e T+1 _i, d { = d T _ mod ^ T) but fa « / T _ 



mod(i,T) 



merely) . 



Generalising (13.71) . the double discrete action has the general form 

^t bl) = k[ d) J2 (Sl du \t) + §W)(t)) + e «iV 4il > - 7V 4 ) 2 . (4.6) 



t=o 



In accordance with (14.41) and the structure of Pdbh I consider the following 
(probably the simplest possible) kinetic term 



Si dU \t) = -a 2 W,i(t+l)-iV4,i(*)) 2 | 



N AA (t) + N 4A (t+V 



2(pN 3 , 2 (t + |) - iV 4jl (t))2 + 2(JV 4 ,i(f) - pN 3:2 (t - \)f 



N 3 , 2 (t-±)+N 4>1 (t) 



a . 



A^ 4j i(t) + aiV 5i2 (t + i 
and the potential one 

S^ dU \t) = -~k 1 2 1 N A>1 (t) 1 s+X n N A , 1 (t)+ 



(4.7) 



(4.i 



where a, a, fc^ 1 ) ^ n > ^i 2 > ^ 22 ( as we ^ as ^ ) are the "double" coupling 
constants and a value of p is taken from the previous section. Consequently, 
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Figure 4.2: -^N^i)-* + A 11 fitted to -d N4 l(i) u. 
trative purposes the term with e is absorbed into A 11 



:) ^ dbl) [iV(t)]; for the illus- 



all elements of Paw depend on the kinetic part of the action S^is while the 



potential part influences merely the diagonal (of the blocks P^, P^y 



The extremum condition (14. 2 p in the expansion (14. ip yields 



Vz: -d 



N p , q (i) 



N, 



(4.9) 



Hence, similarly as in the case of Fig. I3.2[ I may fit first derivative of the 
potential term (14. 8 p to first derivative of the kinetic one ( 14. 7p - for (p, q) = 
(4,1) (Fig. S2]) and (p, q) = (3,2) (Fig. Q]h (The abbreviated notation 
in figures is the same as in the previous chapter.) It is worth to be noted 
that the agreement between them (and hence validity of (I4.2p ) is good in the 
whole relevant region of the universe, almost right up to the stalk and hence 
much better than it was for the single slice structure. As a result, it is well 
justified to examine, making the use of ( 14. 4p (analogically with Section 3.2), 
the substitution of the assumed kinetic and potential term into the second 
term of the expansion ( 14. ip . In particular, the suitability of the kinetic term 
and values of its coupling constants can be directly inferred from 



(-fdbl)v+l 



2e, 



k {d) d 



N 



(4.10) 
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Figure 4.3: \kfN 3 , 2 {i + - A 22 fitted to -d^^S^N®] 

Then, the potential term can be obtained through 

kf ] dl^Sf^Nit)} = (Pj 1 ^ -2s- M d) d^S^imt)], 

*P dl sMl+ ^r ] lN(t)} = (PZh ~ dl s 2{i+l jr ] [N(t)}. (4.11) 

However, fitting LHS's of (I4.10p to their RHS's does not yield completely 
reliable results. Therefore I use another method and determine coupling 
constants one by one. I will not describe that in detail but the scheme is as 
follows. At the beginning I find C x = kf\p + a) 2 from (P]bi)i,i ~ (^dbi)*+i,i- 
Next - knowing the preceding - a from the diagonal (Pdbi)i,v At last, I find 
C*2 = hi a from the superdiagonal (Pdbi)t,t+i- Hence I get a while k[ , a 
can be trivially calculated. Afterwards, the obtained LHS's of f l4.10p may be 
compared to RHS's (Fig. I4.4[ Fig. 14. 5p . with which they should agree. (The 
quite poor result for Fig. 14.41 may indicate the modified >S' j [ dbl ' ) .) Furthermore, 
LHS's of fTCTll with ^ dbl) , ^ dbl) exchanged may be compared to RHS's yet 
with second derivatives of Sp dhl \N(t)] neglected (Fig. l4.6[ Fig. I4.7P . However, 
keeping in mind small magnitude of the latter, such a comparison as well gives 
an insight into the suitability of the kinetic term. At last, one may investigate 
the potential term f 14 . 8 1) using (14. lip and obtaining the concurrent results but, 
for conciseness, I will not show this here. 
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4.3 Reduction to the single slice structure 



Having two different descriptions of the CDT universe, the latter of which 
is intended to be an extension of the former, we would like to be able to go 
in the opposite direction and to reduce the broader framework to the nar- 
rower one. In other words, after "integrating out" (3, 2) simplices we should 
obtain the single slice structure of (4, 1) simplices alone. For the average 
distribution of the number of simplices this is trivial - one just removes the 
second component from the joint distribution of (N^i(i)) and p (N 3 ^{i + \)) ■ 
The situation is the same in the case of the covariance matrix, where simply 
C^i = C^biJ as it was described in the previous section. However, for the 
inverse covariance matrix such a reduction is inevitably not so trivial and 
its existence is by no means obvious. Fortunately, it turns out that for an 
invertible rectangular matrix composed of four square blocks satisfying cer- 
tain basic conditions there are useful algebraic relations between those blocks 
and blocks of the inverse matrix, which we may apply to C^h an d -Pdbi- The 
so-called Boltz-Banachiewicz inversion formula has version I, 

A B 
C D 

(A- 1 + A- l B(D - CA^B^CA- 1 —A~ 1 B{D - CA^B)- 1 } 
\ —(D — CA~ 1 B)~ 1 CA~ 1 (D-CA-'B)- 1 J 1 ^ ] 



and version II, 



A B 
C D 



-l 



(A-BD^C)- 1 -(A - BD^C^BD- 1 

-D^ 1 C{A - BD^C)- 1 D- 1 + D~ 1 C(A - BD^C^BD- 



i > 



(4.13) 



comparison of which yields additionally four algebraic equations. However, I 
will restrict myself to the statement that transforming (I4.13P one can retrieve 
P 4)1 from P dbi : 

Pi,i = {(P^r 1 + CZP S , 2 (OT\ (4-14) 

where the transposition results from the symmetry of Cdbi- 

The next and crucial step is a reduction of the double discrete action 
S^ 1 ^ to the "single" one Sdi S . In the double slice structure I consider the 
kinetic term (14.71) and the potential one (14.81) . the only ones which I have 
discussed in the previous section. Looking at their form one may suspect 
that in the single slice structure they correspond to the (simplest possible) 
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kinetic term ( 13.1 2p and the potential one ( I3.8|) . In order to make the desired 
transition to S^ s I apply a Gaussian approximation to (14.71) i.e. I literally 
integrate out (3,2) simplices, although in a specific way, 

Jdz exp(- (y - x) 2 H- 1 - (z - xfr 1 - (y - zfl' 1 ) = 

constexp(- (y - xf(\H- 1 + r 1 )) (4.15) 

where x = AT 4)1 (t), y = N± 7 i(t + 1), z = pNs^(t + \) and the denominators 
F=-(2fci d) ) :i (iV4,i(t) + iV4 > i(*+l))^= (2^ d) )- 1 (^4,i(*)+^(*+|)) « 
(2A;S d) )- 1 i(l + p(^ r 4,i( t ) + iV 4,i( t + 1 )) (i- e - the approximation is made before 
the integration). For the potential term ( 14. 8 p I assume pN 3t2 (t+^) ~ N^i(t) 

and take the rescaling of implied by (14.15P into account. As a result, the 
"single" coupling constants are expressed through the double ones, 

k 1 = k < f ) (-a + -^-), (4.16) 

k 2 = (-k? + P -r k f) p + ° (4.i7) 

p-a(p + a) 

A = ( - A 11 + p- 1 A 22 ) P + ° - . (4.18) 

The actual investigation produce the following results. Values of the single 
coupling constants obtained directly from the single slice structure (embodied 
in Fig. 13.21 but not Fig. 13. 3[ Fig. I3.4[ where the kinetic term is (I3.15P ) are 
ki = 0.04412 ± 0.00069, k 2 = 17.55 ± 0.87 and A = 0.0543 ± 0.0075 while 
values from the double slice structure, calculated using (I4.16p . (14.17[) . (I4.18p . 
amount to k x = 0.0391 ± 0.0026, k 2 = 14.28 ± 0.84 and A = 0.0530 ± 0.0054. 
The presented uncertainities include only the accuracy of fits and do not 
contain estimates of the systematic errors. 
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Summary 



I have verified the existence of the semiclassical limit of CDT in 3 + 1 dimen- 
sions for spatial hypersurfaces formed by (4,1) simplicial building blocks of 
spacetime. The blob region of the universe generated by the triangulated 
path integral turns out to behave effectively like a discrete mini-superspace 
model. This manifests not only in the average distribution of the number of 
(4, 1) simplices but also in the (quantum) dynamics of the system, governed 
by the discretised mini-superspace action, which is extremised by the aver- 
age distribution. The latter can be investigated using the inverse covariance 
matrix of the number of (4, 1) simplices. Analysis of its exceptionally simple 
structure indicates an intuitive form of the potential part of the action and 
more complicated of the kinetic one. 

Furthermore, I have shown that there is an extension of the above semi- 
classical description - which I called the single slice structure - to the double 
slice structure, including spatial slices of both (4, 1) and (3, 2) simplices. 
After the rescaling of the number of (3, 2) simplices, the joint average distri- 
bution of the number of both types of simplices in the blob again corresponds 
to the discrete mini-superspace solution. Next, the "double" discretised semi- 
classical action, extremised by the average distributions, can be examined us- 
ing the inverse joint covariance matrix, which consists of four simple blocks. 
The only discussed form of the action turns out to work reasonably well. 
Last but not least, the double slice structure can be reduced to the single 
slice structure. 

The next step to make is splitting layers of (3,2) building blocks into 
those of (3, 2) and (2, 3) simplices, which amounts to the triple slice struc- 
ture. I argue for that in Appendix A yet this is a matter of further research, 
the same as e.g. consideration of other forms of the double action. 
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Appendix A 

The triple slice structure 



Below, I present some evidence that the description of the CDT universe from 
Chapters 3, 4 can be extended even more by replacing layers of (3, 2) build- 
ing blocks with the separate treatment of (3, 2) and (2, 3) simplices. Such a 
framework would be the ultimate one because of the eventual inclusion of all 
types of simplices - since (4, 1) building blocks yield the single slice structure 
together with the (1,4) and N^i(i) = Ni^(i), which means (1,4) simplices 
contain no additional information. I call it the triple slice structure, expect- 
ing (3, 2) and (2, 3) simplices to form their own effective spatial slices by 
splitting those of the (3,2), in such a way we will obtain three spatial slices 
in every temporal layer, beginning with spatial tetrahedra which belong to 
the (4,1). Such a split will be in agreement with construction of a causal 
triangulation i.e. the order of gluing of simplicial building blocks presented 
at the beginning of Section 3.1. 

In the first place, let us consider the average temporal distributions of 
the number of (3,2) and (2,3) simplices, using formula ( 12. 2p . Since by defi- 
nition N 3 ^(i) = N 3i2 (i) + N 2j3 (i) individual distributions N 3 k 2 (i), N 2 k 3 (i) are 
reparametrised in the same manner as it was in the case of (3,2) building 
blocks (see Section 4.1). This is justified by the fact that icv's f° r (3,2) 
and (2, 3) simplices differ from those for the (4, 1) to the same degree as the 
(3, 2). Unsurprisingly, I obtain the average distributions (N 3i2 (i)), (N 2)3 (i)) 
having the already known three regions: the stalk (both with 6 simplices at 
least), the blob (with single maxima at % = 39 and i = 40, respectively) and 
the tail. Now, looking at their profiles, I assume (N 3j2 (i)) = (N 3<2 (i + |)) 
and (N 2 ^ 3 (i)) = (N 23 (i + |)). Moreover, fulfilling the condition of invariance 
under time inversion (cf. Sections 2.2, 4.1), I carry out the appropriate proce- 
dures of symmetrisation of the distributions with respect to t = 39 + 1, which 
actually do not make them strictly symmetric for the obvious reason that sets 
{i + |}j and {i + |}j are asymmetric with respect to 39 + |. Nevertheless, 
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(-^3,2 + (-^2,3(^ + |)) are altered only slightly. The invariance under time 
inversion also implies (N 3t2 (i + §)) ~ (N 2>3 (i + §)) ~ \{N 3 ^{i + |)) and it is 
these almost-equalities that explain the obtained shape of the distributions. 
Consequently, in the blob (N 3j2 (i + 3))) (^2,3 (i + §)) should be proportional 
to (N^i{i)) in the interpolated sense (see Section 4.1), with the scaling factor 
2p. They turn out to be indeed. The joint distribution of 2p {N 3t2 (i + |)), 
2p(N 2j3 (i + §)) and (N^ s i(i)) is shown in Fig. \A.1\ along with the average 
fluctuations, rescaled in the case of (3, 2) and (2, 3) simplices. Furthermore, it 
is fitted with the semiclassical solution (I2.5P up to yet smaller deviation than 
in Fig. 14.11 The average total number of simplices (Ns^) ~ (^2,3) ~ 103500. 

Secondly, very briefly, using formula (13 .3p I calculate the covariance ma- 
trices C 3j 2, C 2j 3 and their inverses P 3i2 , -p2,3, which obviously turn out to have 
the same properties as 2 and Ps^, up to a smaller symmetry. Thirdly, also 
very briefly, naturally extending the definition (14.31) I calculate the "triple" 
covariance matrix Ct r i and its inverse P tr \. The triple inverse covariance ma- 
trix P tr i has simple form analogical with P dh \ (cf. (14. 5p ) i.e. it consists of 
nine blocks having band structure. Thus we may expect that it will not be 
difficult to find the corresponding triple discrete action. 

4500 , , , , , , , , , 




t 

Figure A.l: The joint average distribution of (N^i(i)), 2p{N 3 ^{i + |)), 
2p {N 2 ^(i+^)) with the semiclassical solution fitted and the average quantum 
fluctuations denoted by vertical lines (not to be confused with errors!) 
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